function [vN1] = BoundaryCall(S, K, r, T, tau)

vN1 = exp(-r * (tau)) * max((S - K), 0);

